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Abstract. We report the discovery of the optical afterglow (OA) of the long-duration gamma-ray burst GRB 030429, and 
present a comprehensive optical/near-infrared dataset used to probe the matter at different distance scales, i.e. in the burst 
environment, in the host galaxy and in an intervening absorber. A break in the afterglow light curve is seen approximately 1 day 
from the onset of the burst. The light curve displays a significant deviation from a simple broken power-law with a bright 1.5 
magnitude bump with a duration of 2-3 days. The optical/near-infrared spectral energy distribution is best fit with a power-law 
with index j3 = -0.36 + 0. 12 reddened by an SMC-like extinction law with (a modest) A v = 0.34 ± 0.04. In addition, we present 
deep spectroscopic observations obtained with the Very Large Telescope. The redshift measured via metal absorption lines in 
the OA is z = 2.658 ± 0.004. Based on the damped Lya absorption line in the OA spectrum we measure the H i column density 
to be log,/V(Hl) = 21.6 ± 0.2. This confirms the trend that GRBs tend to be located behind very large Hi column densities. 
The resulting dust-to-gas ratio is consistent with that found in the SMC, indicating a low metallicity and/or a low dust-to-metal 
ratio in the burst environment. We find that a neighbouring galaxy, at a separation of only 1'.'2, has z = 0.841 ± 0.001, ruling 
it out as the host of GRB 030429. The small impact parameter of this nearby galaxy, which is responsible for Mg n absorption 
in the OA spectrum, is in contrast to previous identifications of most QSO absorption-selected galaxy counterparts. Finally, we 
demonstrate that the OA was not affected by strong gravitational lensing via the nearby galaxy. 

Key words, gamma rays: bursts - galaxies: distances and redshifts - galaxies: high-redshift - dust, extinction - quasars: 
absorption lines 
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1. Introduction 

Considerable progress has been made in the understanding of 
gamma-ray bursts (GRBs) since BeppoSAX started distribut- 
ing localisations with arcmin precision and van Paradijs et al. 
( 1997) discovered the first optical afterglow (OA). Of the ap- 
proximately 50 OAs that have been detected since 1997, 36 
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have had their redshift spectroscopically determined. As shown 
in Fig. [2 the mean redshift is around 1.4 with 13 GRBs at 
z> 1.5. 

Independently of the brightness of GRB host galaxies, 
OAs can be used as tools to obtain information about the 
gas, metallicity and dust content of the host. In particular, the 
optical/near-infrared (NIR) spectral energy distribution (SED) 
can provide the extinction, while the neutral hydrogen column 
density, A^(H i), can be derived from the damped Lya absorp- 
tion spectroscopy, provided that the burst is distant enough for 
the Lya line to be redshifted into the optical/near-UV domain. 
The detection of damped Lya absorption lines in the spectra of 
several GRB afterglows (Jensen et al. 2001 Fynbo et al l2001bl 
Hjorth et al. 2003a Vreeswijk et al. 2004 1 is consistent with the 
now firmly established link between long-duration GRBs and 
core-collapse of massive stars (e.g. Hjorth et al. 2003b, Stanek 
et al. l2003l . 

Like QSOs, GRB afterglows are distant sources that can be 
used as cosmological beacons to investigate intervening galax- 
ies unrelated to the burst itself. For the current discussion, the 
most important difference between these two astronomical phe- 
nomena is that OAs are ephemeral, vanishing altogether within 
a couple of months. This leaves the line-of-sight clear and with- 
out any interference from a bright object. Another significant 
difference is that afterglows typically have featureless spectra, 
making it easier to identify intervening absorption systems. 

In this paper we present the detection of the OA of 
GRB 030429 and use it to probe the progenitor environment 
(sub-pc scale), the host galaxy properties (kpc scale) and an in- 
tervening Mgn absorber (Gpc scale). The organisation of the 
paper is as follows. The optical, NIR and spectroscopic obser- 
vations are presented in Sect. |2] In Sect. [3] we investigate the 
optical/NIR light curve and the SED. We analyze the spectrum 
and give the redshift of GRB 030429 in Sect.gl In addition we 
fit a H i column density model to a damped Lya- line present in 
the OA spectrum. In Sect. |5] we use the derived properties of 
the OA to compare our results with afterglow models. We also 
discuss the properties of the nearby Mg n absorbing galaxy and 
assess whether it could be responsible for strong gravitational 
lensing of the OA. Finally, the main results are summarised in 
Sect.il 

We adopt a cosmology where the Hubble parameter is 
H = TOkms-'Mpc- 1 , Q m = 0.3 and Q A = 0.7. For these 
parameters, a redshift of 2.66 (0.84) corresponds to a luminos- 
ity distance of 2 1.99 Gpc (5.33 Gpc) and a distance modulus of 
46.7 (43.6). One arcsecond is equivalent to 7.96 (7.63) proper 
kiloparsecs, and the lookback time is 11.1 Gyr (7.0 Gyr). 

2. Observations 

GRB 030429 was detected by the French Gamma Telescope 
(FREGATE), Wide Field X-ray Monitor (WXM) and Soft X- 
ray Camera (SXC) on-board the HETE-2 satellite on 2003 
April 29.446 UT. The burst had a duration of 10.3 s, placing it 
in the "long-duration" burst category. Initially, a 1' radius error 
circle was circulated via the GRB Coordinate Network (GCN) 1 
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Fig. 1. A histogram showing the distribution of GRB redshifts 
known to date (July 2004). The median value of the 36 redshifts 
is indicated by the dashed line, while the mean value is shown 
with the dotted line. 



approximately 2 hours after the burst. A week later the correct 



error radius of 2' was distributed (Doty et al. 2003 ). 

The OA was discovered at the Mt. John 0.6-m telescope 
~3. 5 hours after the burst (Gilmore et al. 2003). It was identi- 
fied outside the initial erroneous SXC 1' error circle after com- 
parison with a DSS-2 red plate. The OA was monitored in the 
optical during the following days with the 1.0-m and 1.55-m 
telescopes at the U.S. Naval Observatory Flagstaff Station 
(NOFS), and with the Very Large Telescope (VLT), using either 
the Antu/FORS 1 or the Yepun/FORS2 combination (Fynbo et 
al. l2003at . In addition, Antu/ISAAC was used to obtain NIR 
images of the OA shortly after the burst. A VLT/FORS2 image 
of the OA and its surrounding field is displayed in Fig.|2] The 
journal of our observations is given in Tabled 

Using Antu we obtained 1 1 x 600 s spectra of the OA start- 
ing on 2003 May 2.214 UT, 2.77 days after the burst (Weidinger 
et al. l2003all2nb3bb . The data were acquired with the FORS1 
instrument in long-slit spectroscopy mode with the G300V 
grism, GG375 order separation filter and a 1'.'3 wide slit. This 
resulted in a wavelength coverage from 3900-8000 A. The slit 
was oriented with a position angle of -71°1 in order to cover 
both the OA and a neighbouring galaxy (lower panel of Fig.|2j. 
The observations were done with the standard resolution colli- 
mator. The individual spectra were combined, yielding a seeing 
of 0'.'8 and a spectral resolution of 10 A FWHM. The projected 
pixel size was 0'.'2 x 2.58 A. 
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Berger & Frail (2003 1 observed the error circle of 
GRB 030429 with the Very Large Array at 8.46 GHz, 2.65, 
5.63 and 16.68 days after the burst. No radio sources were de- 
tected at the position of the OA down to a 3<x limit of 0. 18, 0. 16 
and 0. 10 mJy, respectively. 
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Table 1. A log of the GRB 030429 follow-up imaging obser- 
vations. Upper limits are 2cr in a circular aperture with radius 
1". No correction for Galactic extinction has been applied to 
the photometry. 



At" Telescope/ Magnitude* Seeing Exp. time 



[days] 


Instrument 










[arcsec] 


[s] 


B-band: 
















1.747 


1.0-m NOFS 


22 


24 ± 





30 


2.3 


720 


V-band: 
















0.535 


Yepun/FORS2 


21 


45 ± 





03 


0.7 


5 x 180 


0.764 


Yepun/FORS2 


21 


79 ± 





03 


0.8 


5 x 180 


1.755 


1.0-m NOFS 


21 


89 ± 





32 


2.4 


480 


2.541 


Yepun/FORS2 


23 


18 ± 





11 


0.9 


5 x 180 


6.631 


Yepun/FORS2 




>26 







1.3 


5 x 180 


R-band: 
















0.145 


0.6-m Mt. John 


19 


67 ± 





11 


3.4 


10 x 60 


0.154 


0.6-m Mt. John 


19 


39 ± 





08 


3.3 


10x60 


0.170 


0.6-m Mt. John 


19 


62 ± 





09 


3.3 


10x60 


0.214 


0.6-m Mt. John 


19 


56 ± 





11 


4.8 


10x60 


0.548 


Yepun/FORS2 


20 


86 ± 





04 


0.7 


5 x 180 


0.777 


Yepun/FORS2 


21 


13 ± 





04 


0.8 


5 x 180 


1.761 


1.0-m NOFS 


21 


42 ± 





27 


2.2 


420 


1.862 


1.55-m NOFS 


21 


27 ± 





09 


1.2 


600 


1.884 


1.55-m NOFS 


21 


58 ± 





15 


1.4 


600 


2.553 


Yepun/FORS2 


22 


54 ± 





06 


0.9 


5 x 180 


2.793 


1.55-m NOFS 


22 


55 ± 





09 


1.3 


5 x600 


3.632 


Yepun/FORS2 


23 


71 ± 





12 


0.5 


5 x 180 


6.644 


Yepun/FORS2 


25 


20 ± 





30 


1.2 


5 x 180 


67.641 


Antu/FORSl 




>26 


3 




0.7 


15 x 180 


I-band: 
















0.561 


Yepun/FORS2 


20 


29 ± 





05 


0.7 


5 x 180 


0.790 


Yepun/FORS2 


20 


63 ± 





05 


0.8 


5 x 180 


1.766 


1.0-m NOFS 


20 


51 ± 





37 


2.2 


420 


2.566 


Yepun/FORS2 


21 


72 ± 





09 


0.6 


5 x 180 


6.658 


Yepun/FORS2 


24 


70 ± 





30 


1.1 


5 x 180 


J s -band: 
















0.538 


Antu/ISAAC 


19 


25 ± 





04 


0.6 


20x90 


0.753 


Antu/ISAAC 


19 


51 ± 





04 


0.7 


20x90 


K s -band: 
















0.568 


Antu/ISAAC 


17 


70 ± 





06 


0.6 


30 x 60 


0.783 


Antu/ISAAC 


18 


01 ± 





06 


0.6 


30x60 



" Days after 2003 April 29.446 

* The 0.6-m observations were unfiltered but tied to the i?-band as 
described in the text. 



3. Imaging 

3.1. Photometry 

The data were reduced using standard techniques for de-biasing 
and fiat-fielding. The photometry of the afterglow was car- 
ried out using point-spread function (PSF) fitting photometry 
(Stetson 19871 119941 . For the VLT images, the relative opti- 
cal magnitudes were transformed to the Johnson photometric 
system using observations of Landolt ( 1992 1 standards. For the 
other telescopes we used the calibrated magnitudes of stars in 
the field (Henden[2003 1. We note that the 0.6-m observations 
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Fig. 2. Top: a 180 s fl-band VLT/FORS2 image of the 
GRB 030429 optical afterglow, obtained ~13 hours after the 
burst occurred. The position of the optical afterglow is marked 
with an arrow. The SXC localization is shown with a con- 
tinuous circle of radius 2'. Bottom: a combined 900 s R -band 
VLT/FORS2 image with a seeing of 0'.'5, taken ~3.6 days from 
the onset of the burst. The field of view is smaller than shown 
in the top panel in order to emphasise the surroundings of the 
OA. A galaxy at a distance of 172 from the OA is clearly visible 
and marked with an arrow. 

were unfiltered, but we tied them to the /?-band using Henden 
( 2003 ) standards with similar colours as the OA. For the NIR 
data we used 2MASS stars in the field to transform our obser- 
vations to the standard system. 

As seen in Tabled the afterglow was observed at the first 
two VLT epochs (At * 0.55 days and At as 0.78 days) in all 
five filters (VRIJ S K,) and in VRI during the third VLT epoch 
(Af « 2.55 days). The final detection of the OA is in R and / 
at At x 6.65 days, where Af is the time from the onset of the 
burst. We also have an upper limit in V at the same time, and 
a deep upper limit in R roughly 2 months after the burst. We 
note that there is a galaxy ~172 WNW of the OA, visible in 
the lower panel of Fig. [2] which was initially considered to be 
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Table 2. Photometry of the galaxy J1213. 1-2054.8, 1'.'2 WNW 
of the optical afterglow. No correction for Galactic extinction 
has been applied to the photometry. 



V 


R I 


Js 


K, 


23.70 ±0.18 


22.70 ±0.12 21.29 ±0.09 


20.50 ±0.11 


18.40 ±0.07 



the probable host galaxy (see Sect. 14. 2> . There is no indica- 
tion of a host directly underlying the OA; we can set an upper 
limit of R > 26.3 (2cr) from our last observation. For a future 
host search we include the accurate position of the afterglow 
based on USNO CCD Astrograph Catalog (UCAC2) stars in 
the field: RA(J2000) = 12 A 13 m 07f495(10) and Dec(J2000) = 
-20°54'49'.'54(3) with the error in the last digits indicated in 
parentheses. The accuracy of the afterglow position in RA is 
worse than in declination due to the light contribution from the 
nearby galaxy, which is almost directly west of the afterglow. 

In order to determine the properties of the nearby galaxy 
(hereafter referred to as J1213.1— 2054.8) we used SExtractor 
(Bertin & Arnouts [T9961 to obtain its VRIJ S K S total magni- 
tudes (mag_auto). Given that our first epoch data had a very 
similar seeing in all filters (0'.'6-0'.'7), it was used to extract the 
photometry. The PSF-subtracted images, where the afterglow 
had been removed, were given as input to SExtractor. We used 
the same detection image for the various filters in order to mea- 
sure the flux in an identical aperture. Our results are given in 
Table|2] Note that the galaxy is very red with R - K — 4.3. 

Due to the faintness of the OA at At w 6.65 days, and to 
the fact that the seeing was comparable to the separation be- 
tween the OA and J1213. 1-2054.8, it was not possible to per- 
form simple PSF photometry on the OA. Instead we proceeded 
in the following way. We calculated the OA position relative to 
three nearby stars in the first VLT epoch, where the OA was 
much brighter than J1213. 1-2054.8. We then created a PSF- 
subtracted image (Af « 6.65 days) to remove the majority of 
the J1213. 1-2054.8 light. At the position of the OA we finally 
subtracted the PSF by varying its magnitude until the residuals 
had been adequately removed. 

3.2. Light curve 

Our optical/NIR afterglow light curve is presented in Fig. [3] As 
the R -band was better sampled than the other filters we analyse 
it in detail below. We fit the light curve with a broken power- 
Table 3. The result of fitting the /?-band light curve of 
the optical afterglow with a broken power-law. The differ- 
ent methods are described in Sect. 13.21 We note that x\ Q f = 
X 1 /degree of freedom, is the reduced^ 2 of the fit. 



/Vd.o.f. 


d.o.f. 


at 


«2 


t b [days] 


(i) 6.15 


9 


-0.88 ± 0.03 


-2.87 ±0.25 


2.26 ±0.10 


(ii) 2.33 


6 


-0.99 ± 0.02 


-3.46 ± 0.44 


2.76 ±0.11 


(iii) 1.96 


4 


-0.95 ± 0.03 


-1.72 (fixed) 


0.78 ± 0.01 



t [hours] 

10 100 1000 











- 

□ K s -1 


18 








Js-1 

a 1-1 ; 


20 






A 

...... 


O R 

V+1 
■fcB + 2 


22 




v V 


, "as 




24 












; a, 


= -0.95 ± 0.03 \) 




26 


- a 2 


= -1.72 (fixe 


d) 






: ^ 


= 0.78 ± 0.01 


days 





0.1 1.0 10.0 100.0 

t [days] 



Fig. 3. The optical/NIR light curve of the afterglow based on 
the measurements given in Tabled The dotted line is an ex- 
trapolation of the straight line fit through the early 7?-band data. 
The power-law indices and time of the break refer to the 7?-band 
fit. 



law (with indices a\ and before and after the break, respec- 
tively) in three different ways: (i) including all data points in 
the fit; ( ii) excluding the bump around At « 1.8 days, assuming 
it is short-lived and consisting of three points; (iii) excluding 
the bump, assuming it is long-lived and sampled by six points, 
and fixing at at a value obtained from the SED fitting, where 
at = —p and p is the electron energy power-law index (see 
Sect. 15. 1> . The results of the fitting are listed in Table|3 

Due to the bump in the light curve, ( i) is formally rejected 
with high significance. A better fit is obtained in ( ii) by omitting 
the three points clustered around Af as 1.8 days. On the other 
hand, in this case Aa = ac\ — a% = 2.47 + 0.44, a value too high 
to be consistent with any of the fireball model predictions (see 
discussion in Sect. 15. 1> . The best fit is acquired in (iii), indicat- 
ing there is an increase in the OA flux above the extrapolated 
power-law, lasting for approximately 2-3 days. We note the 
late-time behaviour of the light curve is consistent with the ob- 
servations of Khamitov et al. (2003 ) who find R — 23.5 + 0.5 at 
Af = 3.43 days. 

The other filters display a similar initial decay index, with 

ay ~ ai ~ uks ~ -0.9 and aj s 0.7. In addition, the final 

/-band detection clearly implies that the light curve has steep- 
ened, with a similar late-time decay index as in the /?-band. 
The upper limit in V supports the same conclusion for that fil- 
ter. The bump is also present in the VI filters as indicated by 
the points obtained with the 1 .0-m telescope. The data points 
covering the bump are close enough in time compared to the 
bump time scale in order to compute a reasonable estimate of 
the V — R and R — I colours. These are consistent with the same 
colours calculated from two previous epochs, Af = 0.548 days 
and Af = 0.777 days (see Sect. I3.3i . In conclusion, the light 
curve displays an achromatic behaviour with a bump lasting 
for a couple of days. 
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Restframe wavelength [A] 
7000 5000 4000 3000 2000 1500 



E 



-1.0 






• At = 


0.548 days - 


- 1 .2 




>. 


■ At = 


0.777 days ~_ 


-1.4 










-1.6 










-1.8 




No extinction 

MW 






-2.0 
-2.2 




LMC 
SMC 







14.6 14.8 15.0 15.2 15.4 

LogM [Hz] 

Fig. 4. The spectral energy distribution (SED) of the GRB 
030429 afterglow at At = 0.548 days (filled circles) and 
Af = 0.777 days (filled squares). The data points are based 
on our K S JJRV measurements. The abscissa displays the fre- 
quency/wavelength in the GRB restframe, assuming z = 2.658 
(see Sect. 14. 21 . At both epochs the SED is best fit by an under- 
lying B ^ -0.36 power-law with an SMC extinction law with 
Ay 0.34. 



3.3. Spectral energy distribution of the afterglow 

Our multiband observations of GRB 030429 allowed the 
construction of the SED at two epochs. We interpolated 
the magnitudes to common epochs {At- 0.548 days and 
At = 0.777 days), using the aforementioned power-law indices. 
We note that the flux from the host galaxy has not been sub- 
tracted. However, the host is faint enough (R > 26.3) as seen 
at Af « 67 days, that it should contribute <1% of the flux 
at the epochs we are exploring. The SED was constructed 
as explained in Fynbo et al. (2001ai. The result is shown in 
Fig. 0] where we have corrected the observed data points for 
foreground (Galactic) extinction using the reddening maps of 
Schlegel et al. ( IT9981 . giving E(B - V) = 0.062 at that position 
on the sky. 

In order to quantify the effects of extinction we assume that 
the intrinsic spectrum of the afterglow is a power-law, and fit 
the function F v oc i^xlO - Av to the observed SED, where A v is 
the extragalactic extinction along the line-of-sight to the burst. 
The dependence of A v on v has been parameterised in terms of 
the restframe Ay following the three extinction laws given by 
Pei ( IT9921 for the Milky Way (MW), Large Magellanic Cloud 
(LMC) and Small Magellanic Cloud (SMC). For the assumed 
extinction law, the fit provides B and Ay simultaneously (see 
e.g. Jensen et al. l2U0H Fynbo et al. l2001al Holland et al. l2U0H 
Jakobsson et al. 2003, Hjorth et al. 2003a, Vreeswijk et al. 
2004 1. For comparison purposes we also considered the unex- 
tincted case, a pure power-law spectrum given by F v oc \P. 

The fits are shown in Fig. 0] while the parameters of the 
fits for both epochs are given in Table @] The SMC fit is fully 
consistent with the data from both epochs, while other fits are 



Table 4. The results of fitting different extinction laws to our 
GRB 030429 KJJRV afterglow observations. We note that 
x\ f = X 1 /degree of freedom, is the reduced^ 2 of the fit. 



Fitting function 


2 


d.o.f. 


P 




At = 0.548 days: 










No extinction 


7.45 


3 


-1.22 ±0.04 





MW 


7.84 


2 


-1.37 ±0.07 


<0 


LMC 


8.86 


2 


-0.70 ± 0.24 


0.30 ±0.13 


SMC 


0.14 


2 


-0.44 ±0.17 


0.30 ± 0.06 


At = 0.777 days: 










No extinction 


12.01 


3 


-1.22 ±0.04 





MW 


15.50 


2 


-1.35 ±0.07 


<0 


LMC 


11.08 


2 


-0.31 ±0.24 


0.52 ±0.14 


SMC 


1.17 


2 


-0.27 ±0.17 


0.37 ± 0.06 



formally rejected with high significance. For the redshift of 
GRB 030429 the interstellar extinction bump at 2175 A falls 
into the observed /-band. As expected by the lack of this ab- 
sorption bump in Fig.0J the MW and LMC extinction are com- 
pletely inconsistent with the data. The absence of the 2175 A 
feature, often ascribed to graphite grains and ubiquitous in the 
spectra of sight lines through the Galactic diffuse interstellar 
medium, suggests that the GRB host galaxy has lower metal- 
licity and dust content than the MW and the LMC. In conclu- 
sion, the featureless SMC extinction law provides the best fit at 
both epochs, with an average extinction of Ay = 0.34 + 0.04 
and a spectral index of B = -0.36 + 0.12. This is the most ac- 
curate Ay determination obtained so far for a GRB host galaxy. 
We note that our estimation of Ay is not severely limited by 
the lack of B- and U -band data since the flux observed in those 
bands would be attenuated by the Lya forest. At z = 2.658 (see 
Sect. 14.2b the effective GRB restframe wavelength of these fil- 
ters is indeed located below 1215 A. 

4. Spectroscopy 

4.1. Reduction and spectral extraction 

The data reduction was performed with standard techniques for 
bias and flat-field corrections. The individual spectra were co- 
added, and the resulting combined science frame was sky sub- 
tracted in the following way. In the science frame, regions on 
both sides of the spectra were filtered along the spatial axis with 
a 1 x 13 pixels median filter in order to remove cosmic ray hits. 
These regions were averaged to produce a ID sky spectrum. 
The mean sky spectrum was expanded to a 2D spectrum by du- 
plicating the ID sky spectrum, and this was subtracted from the 
unfiltered science frame. 

The spectra were optimally extracted using the code de- 
scribed in M0ller ( 2000 1. The output of the code is the extracted 
ID spectra and the 2D residual spectra. The spectral PSFs of 
the OA and J1213.1— 2054.8 overlap in our 2D spectrum, so in 
order to separate the two contributions we employed an itera- 
tive procedure: ( i) extract and remove the OA spectrum; ( ii) ex- 
tract the 2D spectrum of the neighbouring galaxy and subtract 
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6000 
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Fig. 5. Left: the normalized OA spectrum 2.77 days from the onset of the burst. We detect five absorption lines in the spectrum. 
The identifications are presented in Table[5] Right: the spectrum of J1213.1— 2054.8. We detect one emission line in the spectrum. 
In both panels the spectra have been rebinned to match the size of a resolution element (10 A), and the lower curve is the Poisson 
error spectrum. Strong residuals from night-sky emission lines are marked with 'n-s'. 



Table 5. Features detected in both the OA and J1213. 1-2054.8 
spectra. The columns show the line identification, the restframe 
wavelength, the observed wavelength, the restframe equivalent 
width (negative for absorption lines) and the inferred redshift 
of the lines, respectively. 



Line 


^rest 


'lobs 


w rest 


z 




[A] 


[A] 


[A] 




Afterglow: 










Lya 






42.1 ±2.9 




Oi/Sin 


1303.3 


4764.8 


2.2 ± 0.4 


2.6559 


Cn 


1334.5 


4885.7 


2.5 ±0.5 


2.6610 


C iv (doublet) 






1.4 ±0.4 




Mg n (doublet) 






3.3 ±0.4 


0.8418 


Jl 2 13. 1-2054.8: 










[On] 


3727.1 


6858.0 


-6.8 ±0.8 


0.8400 



0.5 




-10000 10000 

Relative velocity (km s _1 ) 

Fig. 6. The normalized OA spectrum centred on the Lya ab- 
sorption line. A neutral hydrogen column density fit to the 
damped Lya line is shown with a solid line, while lcr errors 
are shown with dashed lines. 



it from the original spectrum containing both components; (Hi) 
extract the 2D spectrum of the afterglow and subtract it from 
the original spectrum containing both components. After four 
iterations of ( ii) and ( Hi) a stable solution was found. 

The resulting spectra were wavelength-calibrated using the 
dispcor task in IRAF 2 , with an error in the wavelength so- 
lution of the order of 0.3 A (only 3% of our spectral resolu- 
tion). Flux calibration was performed using an observation of 
the standard star LTT7987. Finally, the spectra were rescaled 
to match the V-band magnitude at the same epoch. In Fig.[5]we 
show the spectra of both the OA and J1213. 1-2054.8. 

4.2. Redshift 

We identify three absorption lines at >5cr in the spectrum of 
the OA. We use two of them, namely the narrow lines Oi / 
Sin and Cn, to calculate the redshift of the GRB, placing it 
at z = 2.658 + 0.004. Armed with this redshift we searched 

2 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 



for additional lines at >3cr. At the expected position of the 
C iv doublet we find a likely candidate, although at our spec- 
tral resolution we are unable to resolve it. In the spectrum of 
J1213. 1-2054.8 we detect one emission line, which we tenta- 
tively identify as [On] at a redshift of z — 0.84. Checking this 
redshift against the OA spectrum, we identify an absorption 
line at the location of the Mg n doublet. The redshift of this 
blended Mg n line is based on a minimum x 1 fitting of the ex- 
pected line profile. We conclude that the redshift of the nearby 
galaxy is z — 0.841 + 0.001. The detected lines and their char- 
acteristics are listed in Table|5] These results rule out J1213.1- 
2054.8 as being the GRB 030429 host galaxy. 

4.3. The damped Lya absorption line 

Absorption systems with neutral hydrogen column density 
larger than 2 x 10 20 crrT 2 are referred to as damped Lya ab- 
sorbers (DLAs) as they display strong Lya damping wings even 
in low/medium resolution spectra. These are normally detected 
in QSO spectra but so far four GRB afterglows have shown ev- 
idence for a high column density DLA system (see e.g. Fig. 4 
in Vreeswijk et al. 2004 1. Almost all of the neutral gas in the 
high-z Universe is located in the DLAs and they represent a 
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Table 6. Computation of the closure relation, |ori | + b\[3\ + c, for two afterglow models. A favourable model will have a value 
of zero for the closure relation. The ISM and wind models are for isotropic expansion in a homogeneous and wind-stratified 
environment, respectively. The electron energy power-law index, p, is obtained in two ways, from (3 and ct\ . In addition, Aa is 
calculated from the resulting p in both cases. The numbers in the table are based on case (Hi) in Table|3 



Model 




Vc 


Closure 




p(fi) 


Aa 




p(a, 


) 


Aa- 


ISM 


Vc 


> v a 


0.41 ±0.18 


1-2/3 


= 1.72 + 0.24 


(p + 3)/4 = 1.18 


±0.06 


(4| ai | + 3)/3 = 


2.27 ± 0.04 


1.32 ±0.01 


ISM 


Vc 


< v a 




-2/3 


= 0.72 ± 0.24 












Wind 


Vc 


> v„ 


-0.09 + 0.18 


1-2/3 


= 1.72 + 0.24 


(p+ l)/4 = 0.68 


±0.06 


(4| ai | + l)/3 = 


1.60 ±0.04 


0.65 ±0.01 


Wind 


Vc 


< v a 




-2/i 


= 0.72 ± 0.24 













major fraction of the baryon reservoir out of which stars at 
the present epoch formed (e.g. Storrie-Lombardi & Wolfe 2000 
and references therein). 

In Fig.[6]we plot the normalized OA spectral region around 
Lya. Overplotted is a fit to the strong Lya absorption line 
yielding logA^(Hi) = 21.6 + 0.2, within the range of previ- 
ously known GRB-DLAs (21.2 < logN(Hi) < 21.9). The 
redshift deduced in the previous subsection, z = 2.658, was 
imposed on the fit. Assuming Ry = 2.91 and using Ay = 
0.34 + 0.04 from the best fit extinction law (SMC), the ratio 
between the neutral hydrogen column density and the redden- 
ing in the GRB 030429 host is N(H i)/E(B - V) = (34 + 16) x 
10 21 crrT 2 mag -1 . We note that this is a lower limit as J 12 13.1- 
2054.8 could contribute to the measured Ay. This value is con- 
sistent with that of the SMC, (44 + 7) x 10 21 cirT 2 mag -1 , and 
similar to those found for other GRB absorbers (see Table 3 in 
Hjorth etal. l2003al . 

There is an indication of Lya emission in the centre of the 
trough, but it is not statistically significant (< 2cr). The flux of 
the line, if real, is approximately 2.5 x 10~ 18 ergs -1 cm -2 . In 
our assumed cosmology this corresponds to a luminosity of 
1.5 x 10 41 ergs -1 . We can use this result to estimate the star- 
formation rate (SFR), assuming the conversion from Lya lumi- 
nosity to the SFR of 10 42 erg s -1 = 1 M yr -1 (Kennicutt fT998l 
Cowie & Hu fT998l . The Lya SFR in the host of GRB 030429 is 
thus ~0.15M o yr -1 . If confirmed this will strengthen the con- 
clusion of Fynbo et al. ( 2003b I that Lya emission is much more 
frequent among GRB host galaxies than among the Lyman- 
Break galaxies at similar redshifts and that there could be 
a low metallicity preference for GRBs (see also Woosley & 
MacFadven lT999l 



5. Discussion 

5.1. Comparison with afterglow models 

The nature of the ambient medium in which the GRB orig- 
inated can be probed with the parameters <x\, ao and B. We 
consider two afterglow models: an isotropic expansion into a 
homogeneous medium (Sari et al. 1998 1, and an isotropic ex- 
pansion into a wind-stratified medium (Chevalier & Li 1999 1. 
We use the closure relation introduced by Price et al. ( 2002 1 
in order to differentiate between the models. In our notation 
|ai| + b\B\ + c = 0, where the values of b and c depend on the 



location of the cooling frequency, v c , relative to the optical/NIR 
bands, v D , at the epoch of the observations. 

In Table [6] we have used the average value of the spectral 
index,/? = -0.36 + 0.12, to derive the electron energy power- 
law index, resulting in p — -2/3 = 0.72 + 0.24 if v c < v Q and 
p = 1 — 2/8 = 1.72 ± 0.24 if v c > v Q . The former value is 
not considered relevant in afterglow models, and the accurate 
analytical spectral and temporal indices have not been derived 
in the literature for this case. 

The ISM model is not favoured by the closure relation. We 
have used the results from case ( Hi) (see Sect. 13. 2\ . but note that 
( i) and ( ii) give a similar result as a\ does not vary significantly 
between all three cases. 

Each spectral or temporal power-law index relates to a 
certain value of the electron energy power-law index, p, and 
the correct model should result in a similar p for all indices. 
Assuming that the fireball model is valid, this simple fact ex- 
cludes (i) and (ii). In these cases p(J3) + p(ai) + p{a\), and the 
predicted Aa is in addition approximately twice as small as the 
observed one. 

In conclusion, only the wind (v c > v Q ) model produces 
a closure relation consistent with zero. We have fixed ai = 
-p(J3), case (Hi), which results in an observed Aa = 0.77+0.03. 
This is marginally consistent with the predicted values listed in 
Table|6] 

The value of p as estimated from /3, p(J3) = 1 .72 + 0.24, 
is marginally consistent with p - 2. We derived the closure 
relations for the flat spectrum case, 1 < p < 2 (Dai & Cheng 
2001 Bhattacharva I2001i . and found that although formally 
consistent with the observationally determined value of p(J3), 
none of the closures were consistent with zero. In addition, the 
predicted Aa values were quite far from that observed. This 
suggest that the p > 2 case is a better representation of the data, 
although not perfect. This may be due to insufficient sampling 
of the light cure and its overall structure (e.g. bumps). 

Accepting the above scenario, it is clear from Fig. [3] that 
there is a bright long-lived achromatic bump above the late- 
time light curve fit. Due to the sparse data sampling we can- 
not constrain the properties of this bump in detail, but it lasts 
for 2-3 days and seems to peak at 1.5 magnitudes above the 
fit. Similar undulations have been observed previously in light 
curves associated with GRBs: GRB 970508 (e.g. Galama et 
al. H9951 Sokolov et al. [rWSl, GRB 000301 C (e.g . Sagar et 
al. l2000l Garnavich et al. l2000l Masetti et al. 2000 Berger et 
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Fig. 7. Top: the reduced^ 2 as a function of the photometric red- 
shift of the nearby galaxy, obtained from hyperz. The best fit is 
obtained for z = 0.82. The spectroscopic redshift of 0.841 is 
indicated by the vertical dashed line. Bottom: the best fit tem- 
plate galaxy spectrum along with the observed VRIJ S K S SED 
based on the measurements given in Table [2] (filled circles). 



al. EOOOl Jensen et al. EOOTl Rhoads & Fruchter l2"U0Tl Dai & 
Lu l200Tl Panaitescu l200Tl Gaudi et al. l200"!l . GRB 000911 
(Ramirez-Ruiz et al. l2002l . GRB 011211 (e.g. Holland et al. 



2002 Jakobsson et al. 2004i, GRB 021004 (e.g. Lazzati et al 



2002 Bersier et al. 2003 Heyl & Perna 2003 Holland et al. 
I5U0H Mirabal et al. 1277011 Fox et al. lSUOll and GRB 030329 
(e.g. Lipkin et al. 2004 and references therein). 

Deviations from the common smooth power-law decay 
have been explained by a variable external density (e.g. Wang 
& Loeb 2000 1, refreshed shocks from the inner engine (e.g. 
Rees & Meszaros 119981 . a non-uniform jet structure (e.g. 
Meszaros et al. 1998 1 or a microlensing event (e.g. Garnavich 
et al. I2000> . Simultaneous multi-wavelength observations can 
help to determine which of the above scenarios is responsible 
for OA light curve variations (e.g. Nakar & Piran 2003 , Nakar 
et al. 2003, Jakobsson et al. 2004 1. Armed only with a sparse 
optical/NIR data set for GRB 030429, we are unable to distin- 
guish between these possibilities. 



5.2. The intervening absorption galaxy 

In order to gain more information about J1213.1— 2054.8 we 
used the photometric redshift code hyperz (Bolzonella et al. 
2000 1. In the fitting process the programme uses standard^ 2 
minimisation, i.e. computing and minimising the deviations be- 
tween the observed SED (see Table |2j of an object and tem- 
plate SEDs compared within the same photometric system. For 
J1213.1— 2054.8 the best fit is obtained for a photometric red- 
shift of z = 0.82+0.02, in excellent agreement with the spectro- 
scopic redshift of 0.841 that we previously determined (upper 
panel of Fig. 0). 

The best fit template spectrum is displayed in Fig. 
(lower panel). It corresponds to a fairly young stellar popu- 
lation, age between 0.5-1.0Gyr, with a significant extinction 
of Ay = 0.8. In addition, in J1213. 1-2054.8 we have detected 
[O n] in emission (Tabled. All evidence is thus consistent with 
a star-forming galaxy. We calculate its absolute restframe R- 
band magnitude by transforming the observed /-band magni- 
tude (Table |3} to AB magnitudes (Fukugita et al. [1995 1, ap- 
plying the distance modulus and the 2.51og(l + z) term (see 
e.g. Appendix A in van Dokkum & Franx 1996). This results 
in Mr(AB) = -21.6, a value located on the bright end of the 
luminosity function. 

The galaxy is located 1'.'2 from the line-of-sight to the OA, 
which corresponds to an impact parameter of D — 9 kpc at 
z = 0.841 in our adopted cosmology. D is frequently inter- 
preted as the lower limit of the radii of suspected counterparts 
of QSO absorption-selected galaxies. Guillemin & Bergeron 
(1997) reported that galaxy counterparts of Mgn absorbers 
with equivalent widths larger than 0.6 A and 0.7 < z < 1.3, 
all had D > 50 kpc. That result was consistent with the conclu- 
sions of Bergeron & Boisse ( 1991 ), who examined a sample of 
z ~ 0.5 Mgn absorbers (see also SteidelQ995j). 



In every case, the possibility exists that the real absorber is 
very faint and/or hidden inside the glare of the QSO, thus over- 
estimating the suspected counterpart radius. Indeed, more than 
half of present-day galaxies reside in groups (Eke et al. 2004), 
lending support to this hypothesis. The fact that many OAs are 
high-z (Fig.Q and transient phenomena has made it possible 
to search for faint/nearby absorbers. The OAs fade away com- 
pletely within a few months, so it is only a matter of obtaining 
long exposures in order to get deep limits. GRB 030429 is the 
second burst in which a Mg n absorber with a small D has been 
identified with imaging and spectroscopy. The GRB 020405 
sight line also revealed a spectroscopically confirmed interven- 
ing Mgn absorber with an impact parameter of only 13 kpc 
(Masetti et al. l2"0"0"3l . 

In Table we list all GRBs that have an Mgn absorption 
system along the line of sight in addition to that originating 
from the host galaxy. We have omitted bursts whose spectra 
display an Mg n absorption that most likely arises in a system 
physically interacting with the host (GRB 000926: Castro et 
Barth et al. 2003 ). Of the seven bursts, 



al. 120031 GRB 020813: 
two (GRB 020405 and GRB 030429) have a spectroscopically 
confirmed counterpart with a very small impact parameter of 
D ~ 10 kpc. The host of GRB 991216 is located close to at 
least three nearby galaxies which are likely candidates for the 
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Table 7. GRBs whose spectra display an Mg n absorption sys- 
tem besides the one produced in the host galaxy. Here, D is the 
impact parameter, the projected distance between the OA and 
the candidate Mg n absorber. The impact parameter is given in 
kpc where the absorber is spectroscopically confirmed. 



GRB 


ZGRB 


<-abs 


D 


970508 


0.835 


0.768 


~5'.'5 


991216 


1.022 


0.770, 0.803 


~0'.'5, ~2'.'0 


010222 


1.477 


1.156, 0.927 


~2'.'0, ~4'.'0 


020405 


0.691 


0.472 


13 kpc (2'.'0) 


021004 


2.335 


1.381, 1.604 


-16" 


030226 


1.986 


1.043 


~2'.'5 


030429 


2.658 


0.841 


9 kpc (1'.'2) 



two Mgn absorption systems (Vreeswijk 2002). Should that 
be the case, their impact parameters would correspond to only 
4-15 kpc. The situation is similar for GRB 010222, where at 
least four nearby galaxies surround its host (e.g. Mirabal et al. 
120021 Frail et al. l2002l Galama et al. 2003). Two galaxies are 
located 2Y5 away from GRB 030226, corresponding to 20 kpc 
at z = 1.043 (e.g. Price et al. 120031 Klose et al. 2004 1. A pos- 
sible galaxy counterpart to the Mg 11 system in GRB 970508 
is located 5'.'5 away from the host, corresponding to 40 kpc at 
z = 0.786 (e.g. Metzger et al. [T9971 Djorgovski et al. [79971 
Pian et al. lT39lfl Sokolov et al. fT^l Fructher et al. lIUUOl For 
GRB 021004 a possible galaxy candidate at an impact param- 
eter of 16", corresponding to 135 kpc at z = 1.381, has been 
observed in narrow-band imaging (Vreeswijk et al. 2003 1. The 
authors point out that a deeper study should be performed on 
the GRB 021004 field before a faint absorber at a low impact 
parameter can be ruled out. 

A large fraction of possible OAs appears to have a faint ab- 
sorber at small D. These results imply that some previous QSO 
absorption-selected galaxy counterpart identifications may be 
incorrect, with the real absorber being a nearby faint galaxy. 
This suggestion is supported by the observations of Falomo et 
al. ([19971 and Watson et al. (2004i, who find a likely counter- 
part to an Mg 11 absorption system in a BL Lac object at only 
D = 11 kpc. In addition, Ellison et al. (2004i analysed 27 in- 
tervening absorption systems in the spectra of a triply imaged 
QSO. They find that the most likely absorber coherence scale 
is 3 kpc, more than one order of magnitude smaller than the di- 
mensions of Mg 11 absorbers deduced from the impact parame- 
ters. 

5.3. Strong gravitational tensing? 

With an angular separation of only 1'.'2, J1213. 1-2054.8 could 
affect the appearance of the OA via strong gravitational lens- 
ing (e.g. Narayan & B artelmann 119991 Grossman & Nowak 
1994 1. For that to occur the Einstein radius must be at least 
1'.'2, which corresponds to a lens mass of>5.4xl0 n M o . Using 
M S (AB) = -21.6 gives a mass-to-light ratio of >15 T (where 
T Q = Mq/Lq). At z — 0.841 we are only probing the innermost 
9 kpc of the galaxy. The aforementioned mass is roughly six 
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times larger than for Milky Way-like galaxies within a simi- 
lar radius. In addition, the critical surface mass density for this 
configuration is 0.43 gem -2 , a factor of three to four times that 
in normal galaxy lenses (see discussion regarding GRB 990123 
in Andersen et al. U999> . 

We also utilised the lensmodel package (Keeton 2001 1 to 
model J 12 13. 1-2054.8 as a singular isothermal ellipsoid model 
with a flat rotation curve. Using our /?-band image from At = 
67.6 days we constrained its ellipticity, e x 0.16, and posi- 
tion angle, <p ~ -49°. Knowing the distance between the lens 
and the OA, we varied the lens mass until the f indsrc output 
resulted in a multiple image system. This occurred when the 
mass inside 9 kpc reached >3 x 10 11 M , a value comparable 
to the one obtained above. In conclusion, the three independent 
methods discussed here virtually rule out that GRB 030429 was 
multiply imaged (with a time delay of roughly a month) by the 
nearby galaxy. 



6. Conclusions 

GRB 030429 occurred in a faint (R > 26.3) galaxy at a redshift 
ofz = 2.658+0.004. The derived neutral hydrogen column den- 
sity is logAf(Hi) = 21.6 + 0.2, while the restframe reddening, 
obtained from the optical/NIR SED, is accurately measured to 
be A v = 0.34 + 0.04. We conclude that the high value of the 
ratio between column density and optical extinction is readily 
explained by a metal-poor environment, similar to that of the 
SMC. This is fully consistent with the favoured SMC-like ex- 
tinction law obtained from the SED fit. 

A nearby galaxy with a separation of only 172 is ruled out 
as the host, with z = 0.841 + 0.001. This is the second time 
a galaxy adjacent to a GRB OA is found to be responsible for 
a Mg 11 absorption system in the OA spectrum. In both cases 
the impact parameter was small, D ~ 10 kpc. At least three ad- 
ditional bursts have galaxy counterparts with similar small im- 
pact parameters that presumably account for their Mg 11 absorp- 
tion systems. Thus, at least five out of seven GRB OAs display- 
ing Mg 11 in absorption (see Table[7), have a nearby galaxy with 
D ~ 10-20 kpc. This strongly indicates that previous identifi- 
cations of many QSO absorption-selected galaxy counterparts 
are possibly incorrect. The remaining OA galaxy counterparts 
should be spectroscopically confirmed before a firmer conclu- 
sion can be drawn. 

Within the framework of the afterglow synchrotron model 
(Sari et al.[T998 1, the GRB 030429 early-time decay index, a u 
and spectral index, (3, indicate that a jet (v c > v D ) expanding 
into a wind medium provides the most likely scenario for the 
data. The late-time decay index, at, has to produce an electron 
energy power-law index, p, equivalent to that calculated from 
/}. This suggests that a bright long-lived bump is present in the 
GRB 030429 light curve. We emphasise that densely sampled 
GRB light curves are of the utmost importance in order to ob- 
serve these variations; without the At » 1.8 days GRB 030429 
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observations it would have been difficult to detect the flux in- 
crement above the interpolated power-law. 

Note added in manuscript. — After the submission of this pa- 
per, GRB 030429 was classified 3 as an X-ray flash (XRF). 
Together with the unambiguous OA redshift determination (see 
Sect. 14.21 . this makes it the first XRF with a spectroscopically 
confirmed absorption redshift. We note that within the uncer- 
tainties the burst is consistent with being an X-ray rich GRB. 
The most probable reason for this burst being classified as an 
XRF is a combination of its fairly low intrinsic £p e ak, and the 
relatively high redshift of the burst (resulting in an observed 
£peak ~ 35keV). 
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